arXiv:1502.00663vl [math.AP] 30 Jan 2015 


STABILITY OF OBSERVATIONS OF PARTIAL DIFFERENTIAL EQUATIONS 
UNDER UNCERTAIN PERTURBATIONS 

MARTIN LAZAR 


Abstract. We analyse stability of observability estimates for solutions to wave and Scrodinger equations 
subjected to additive perturbations. The paper generalises the recent averaged observability/control result 
by allowing for systems consisting of operators of different types. The method also applies to the simultaneous 
observability problem by which one tries to estimate the energy of each component of a system under 
consideration. The analysis relies on microlocal defect tools; in particular on standard H-measures, when 
the main dynamic of the system is governed by the wave operator, while parabolic H-measures are explored 
in the case of the Schrodinger operator. 


1. Introduction 

A notion of averaged control has been recently introduced in m Ho] , both for parameter dependent 
ODEs, as well as for systems of PDEs with variable coefficients. Its goal is to control the average (or more 
generally a suitable linear combination) of system components by a single control. The problem is relevant 
in practice where the control has to be chosen independently of the coefficient value. 

The notion is equivalent to the averaged observability, by which the energy of the system is recorded by 
observing the average of solutions on a suitable subdomain. 

In the paper we investigate a more general problem based on a system whose first component represents 
its main dynamic, while the other ones correspond to perturbations. Assuming that the main component is 
observable, we explore conditions by which that property remains stable under additive perturbations. 

In general, the operators entering the system are not assumed to be of the same type. In a special case 
of a system consisting of a same type operators, the result corresponds to the averaged control of the system, 
thus incorporating results obtained in [10] . 

The methods are applied to the simultaneous observability problem as well, by which one tries to estimate 
energy of all system components by observing their average. The corresponds dual problem consists of 
controlling each individual component of the adjoint system by means of a same control. 

The study of the problem explores microlocal analysis tools, in particular H-measures and their variants. 
H-measures, introduced independently in [71[T3], are kind of defect tools, measuring deflection of the weak 
from strong convergence of L^ sequences. Since their introduction, they have been successfully applied 
in many mathematical fields - let us just mention generalisation of compensated compactness results to 
equations with variable coefficients mm and applications in the control theory maun]. Most of these 
applications apply the so called localisation principle providing constraints on the support of H-measures 
(e.g. [2]), and the proofs of this paper rely on it as well. 

The paper is organised as follows. In the next section we provide an averaged observability result for a 
system whose main dynamic is governed by the wave operator. The finite system is analysed first, followed 
by generalisations to an infinite discrete setting. Application of the approach to simultaneous observability 
is provided in the subsection 2.3. The third section is devoted to observation of the Schrodinger equation 
under perturbations determined either by a hyperbolic or by a parabolic type operator. In the latter case, 
parabolic H-measures (generalisation of original ones to a parabolic setting) have to be explored. The paper 
is closed with concluding remarks, and by pointing toward some open and related problems. 
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2. Observation of the wave equation under uncertain perturbations 


2.1. Averaged observability. We analyse the problem of recovering the energy of the wave equation by 
observing an additive perturbation of the solution. The perturbation is determined by a differential operator 
P 2 , in general different from the wave one. 

More precisely, we consider the following system of equations: 


Pi Ml 


div (Ai(t,x)VMi) = 0, 

(t,x) G R"*" X H 


P 2 U 2 = 0, 

(t, x) G R^ X H 


Ml = 0, 

(t, x) G R^ X dll 

(2.1) 


Mi(0,-) = /3o e 

dtMi(0, •) = /3i e , 


where is an open, bounded set in R'^, Ai is a bounded, positive definite matrix field, while P 2 is some, 
almost arbitrary, differential operator (precise conditions on it will be given below). In the sequel we shall 
also use the notation Ai = —div (AiV) for the elliptic part of Pi. 

For the moment, we specify neither initial nor boundary conditions for the second equation, we just 
assume that corresponding problem is well posed and that it admits an solution. For the coefficients of 
both the operators we assume that are merely bounded and continuous. 


Proposition 2.1. Suppose that there is a constant C, time T and an open subdomain uj such that for any 
choice of initial conditions /3q,/3i the solution ui of jni) satisfies 

Pi(0) := \\Ml 2 + m\l-i <C r [ l^ipdxdt. (2.2) 

Jo j LO 

In addition, we assume that characteristic sets {pi(t, x, t, ^) = 0},i = 1,2 have no intersection for (t,x) € 
(0, T) X uj, (r, ^) G where pi stands for the principal symbol of the operator Pi. 

Then for any 61,62 & 'R.,9i ^ 0 there exists a constant such that the observability inequality 

Pi(0)<C,, ^ y |diMi + 02U2pdxdt+||/3o|lH-i+ II^i|Ih-^ (2-3) 

holds for any pair of solutions (mi,M 2 ) to (12.11) . 


Proof: Of course, the case 62 = 0 holds trivially, and for simplicity is excluded from the further analysis. 
We argue by contradiction. Assuming the contrary, there exists a sequence of solutions u" ,it 2 such that 

Pi"(O)>n^^y|0i< + 02i^^pdxdt+||^o”llH-i + ll/3rilH^^ ■ (2.4) 

As the considered problem is linear, without loosing generality we can assume that if"(0) = 1. Thus (12.41) 
implies that ||,5 oIIh-i + ll/^rilH-^ 0, resulting in the weak convergence (/Jg,/?") —^ (0,0) in L^(r2) x 
H“^(0). Therefore the solutions (u") converge weakly to zero in L^(n x (0,T)) as well. In order to obtain 
a contradiction, we have to show that the last convergence is strong, at least on the observability region. 

From the contradictory assumption (12.41) we have that the H-measure u associated to a subsequence of 
(01 u" + 62 U 2 ) vanishes on {0,T) x uj. Furthermore, it is of the form 

u = diHi + 021*2 + 0i0225i/ri2, 

where on the right hand side the elements of the matrix measure associated to the vector subsequence of 
(mi,M 2 ) are listed, with fii 2 denoting the off-diagonal element. Note that (rt^) is bounded in L^((0,r) x w), 
since that is the case for (m") (by boundedness of initial data), and for the linear combination (0 iMi -I- 62 U 2 ) 
(by contradictory assumption (12.41) 1. which enables one to associate an H-measure to it. 

According to the localisation property for H-measures, each ptj is supported within the corresponding 
characteristic set {pi(t,x, r, ^) = 0}, z = 1,2, which, by assumption, are disjoint on the observability region. 
On the other hand, from the very definition of matrix H-measures it follows that off-diagonal entries are 
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dominated by the corresponding diagonal elements. More precisely, it holds that supp ^12 ^ supp fl 
supp ^ 2 , implying that /ii 2 = 0 on the observability region. 

Thus we get that 

u = + 02^2 = 0 on (0,r)xa;. 

As fii and ^2 are positive measures and 9i > 0, it follows that vanishes on (0, T) x w as well. Thus 
we get strong convergence of (u") in L^((0,T) x ui), which together with the assumption of the constant, 
non-zero initial energy contradicts the observability estimate (12.21) . □ 

Remark 1. The last result provides surprising stability of the observability estimate (12.21) under uncertain 
perturbations, up to compact reminders. Essentially, the only requirement for the perturbation is separation of 
the characteristic sets. This implies that the wave component can be observed robustly when adding unknown 
perturbations, up to a finite number of low frequencies. 

In the next step we would like to obtain the strong observability inequality for initial energy i?i(0) by 
removing compact terms in (12.31) . To this effect, we have to specify some additional constraints on the 
problem for the perturbation U 2 . 

We take P 2 to be an evolution operator of the form 

P 2 = {dt)'^+ C 2 {^)Ai, kGN, (2.5) 

where Ai is an elliptic part of the wave operator Pi, while C 2 is a bounded and continuous function. 

Theorem 2.2. As above, we assume that the coefficients of the operator Pi are bounded and continuous, 
and that the corresponding solution ui satisfies the observability inequality (12.21) . 

In addition we assume that the perturbation operator P 2 is of the form (ES. In the case k = 2 (i.e. P 2 
being a wave operator) the separation of coefficients C 2 (x) — 1 0 is supposed on tv. 

For the initial values of solutions Ui,i = 1,2 we supposed these are related by a linear operator such that 
whenever ((0iWi(O) -I- 02 U 2 {O))\ui = 0) then (ui(0)|^ = M2(0)||,^ = O), and the analogous implication holds for 
the initial first order time derivatives. 

Then there is a positive constant such that the strong observability inequality holds: 

Ei{ 0 )<Csi ( f \ 6 iui + 92U2\^ dxdt. ( 2 . 6 ) 

Jo J uj 

Remark 2. Note that the above assumptions directly imply that characteristic sets of Pi and P 2 are disjoint. 
Indeed, for P 2 being an evolution operator of order k ^ 2, its principal symbol equals 0 only in poles ^ = 0 
(case k = 1), or on the equator r = 0 (case k > 2 ) of the unit sphere in the dual space, where pi = 
— Ai(x)^^ differs from zero. 

In the case k = 2 separation of the characteristic sets is provided by the assumption C 2 (x) 1 on u. 

Proof: As in the proof of Proposition 12.11 let us suppose the contrary. Then there exists a sequence of 
solutions u”,it 2 to (12.11) such that A"(0) = 1 and 

[ [ \0iu^+ 92ulf\^dxdt—^0. (2.7) 

Jo J UJ 

Thus the corresponding weak limits satisfy both the equation PiUi = 0, as well as the relation 

0lUi -f 02U2 = 0 (2.8) 

on the observability region. As functions Ui,dtUi are continuous with respect to time, the bounds on initial 
data imply /3o = /3i = 0 on w. 

Assumptions on the operators Pi,i = 1,2 ensure that corresponding characteristic sets do not intersect. 
By applying localisation property of H-measures as in the proof of Proposition 12 .1 1 we get that u" converges 
to ui strongly on (0,r) x u. 

It remains to show that the limit ui vanishes on the observability region, which, together with the 
assumption of the constant non-zero initial energy, will contradict the observability assumption (12.21) . 

We split the rest of the proof into several cases. 
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a) (k = 2) Due to the relation (12.81) it follows that 

(c 2 - = 0, (t,x) e (0,r) X w. 

As |c 2 — 1| > 0, and the initial data are 0 on uj, it implies tti = 0 on the observability region. 

b) (k = 1) Relation (j2.8|) implies 

dtUl = C2(x)9ttMi, 

which together with ui(0, •) = dtUi{0, •) = 0 on w provides the claim. 

c) (k > 2) Similar as above we obtain 

dt~‘^{dttui) = C2(x)5ttui . 

As •) = —yliUi(0, •) = 0 on w, and similarly for the higher order derivatives, the claim 

follows. 

□ 


Remark 3. Several remarks are in order. 

• The observability assumption (12.2|) on a solution of the wave equation is equivalent to the Geometric 
Control Condition (GCC, [3]j, stating that projection of each bicharacteristic ray on a physical space 
has to enter the observability region in a finite time. 

• The last theorem also holds if, instead of initial data of two components being linked by an operator, 

we assume a cone condition \\u 2 < c||wi( 0 )||l 2 (,^) (||m 2 ( 0 )||l 2 (^) > c||mi( 0 )||l 2 (,^)), with a 

constant c < 0i/02 (c > The latter condition is stable under passing to a limit, and also 

ensures the implication ((0iUi(O) + 92U2{0))\ui = 0) => (ui(0)|(^ = U2(0)||;^ = O), which suffices for 
the proof. 

• The result (12:61) can be generalise to a more general perturbation operator P 2 by assuming that 
coefficients of both the operators are analytic. In that case the separation of characteristic sets 
implies the separation of corresponding analytic wave front sets m Theorem 9.5.1]/ Together with 
(12.81) it provides that ui is analytic on the observability region. Constraints on initial data and finite 
velocity of propagation imply ui = Q on an open set near t = 0 , and as the solution is analytic it 
vanishes on the whole observability region which contradicts the observability assumption (12.21) . 

• If P 2 is a wave operator the strong observability inequality (ESI) is equivalent to the controllability of 
a suitable linear combination (determined by the operator linking the initial data) of solutions to the 
adjoint system under a single control (cf. [10] ). 

Meanwhile, the weak observability result (|2.3I) in that case corresponds to the average controllability 
of the adjoint system up to a finite number of low frequencies. 

• The last theorem generalises the results of [10] by allowing for a general evolution operator P 2 which 
does not have to be the wave one. 

In addition, it allows for an arbitrary linear combination of system components, while in [lOj 
just their (weighted) average is explored. Specially, if the difference ui — U 2 is considered, the result 
corresponds to the synchronisation problem (e.g. [12] ) in which all the components are driven to the 
same state by applying the null controllability of their differences. 

Furthermore, unlike in [lOj , the proof of the relaxed observability inequality (|2.3p does not rely on 
the propagation property of H-measures, which allows for system’s coefficients to be merely continu¬ 
ous. On the other hand, such approach avoid technical issues related to the reflection of H-measures 
on the domain boundary. 

• The theorem to some extent also generalises the results of [m in which a similar result is provided 

for the system consisting of a wave and a heat operator with constant coefficients (or more 

generally with a common elliptic part). 

However, although allowing for a more general perturbation operator, it requires initial data of 
two components of the system to be related, while in [E] no assumptions on initial data for the 
second component is assumed. 

• The constant can be taken to be uniform for 61 > 0* and 0* > 0. 
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• The weak observability result (lOl) is easily generalised to a system with a finite number of compo¬ 
nents, under assumption that the characteristic set of the leading operator Pi is separated from the 
characteristic sets of all the other operators, while the latter ones can be arbitrary related. 

However, the generalisation of the result to an infinite dimensional setting is not straightforward. 
It requires study of the localisation property for H-measures determined by a sequence of function 
series, and is the subject of the next subsection. 

On the other side, the generalisation of the strong estimate result (1^ to a system consisting of 
more than two components has still not been obtained, and is a subject of the current investigations. 

• The result (12311 also holds if the observability region is not of a cylindrical type, but a more general 
set satisfying the GCC for the first component of the system. Such generalisation corresponds to a 
moving control (cf. [11] ). However, possible derivation of the corresponding stronger observability 
result remains open, due to the last part of the proof of Theorem \2.‘A in which the special (cylindrical) 
shape of the observability region is used. 

• Note that the observability result (1^ is weaker than the one required in the simultaneous control, 
where one has to estimate initial energy of all components entering the system. Of course, the 
assumptions in the latter case are stronger, as one has to assume that the observability set satisfies 
the GCC for the second component as well. 


2.2. Infinite discrete setting. In this subsection we want to analyse the stability of the observability esti¬ 
mates for the wave equation when the perturbation is given as a superposition of infinitely many components, 
each determined by a differential operator Pi, which, in general, does not have to be the wave one. Thus the 
system of interest reads as: 

Piui = dttui — div(Ai(t,x)VMi) = 0, (ijx) s R+ x O 

PiUi = 0, {t, x) £ R'*' X fl,i > 2 

ui=0, (t,x)£R+x5n (2.9) 

ui(0,-) =/3o eL2(r!) 

dtui{0,-) = /3i G , 


with the same assumptions on the domain and the operator Pi being assumed for the system iH). For 
the other equations, neither initial nor boundary conditions are specified. For the moment, we just assume 
the corresponding coefficients are bounded and continuous, and the problems are well defined with solutions 
in L^(n). 

In this setting, the same microlocal analysis tool as in the finite case, in particular the localisation 
property of H-measures, is applied in the study of the stability of the observability estimates. However, as 
perturbations are determined by a superposition of infinitely many solutions, this requires analysis of the 
mentioned property for a sequence of function series. 

Namely, it is well known that an H-measure associated to a linear combination of two sequences is 
supported within the union of supports of measures determined by each component, and the same property 
holds for any finite linear combination. However, in general it fails when considering superposition of infinite 
many sequences, as shown by the next example. 


Example 1. Let (u") and (/”) be L^(R'^) sequences, whose corresponding H-measures pu and pf have 
disjoint supports, and let {Of) be a sequence of nonnegative numbers summing in 1. 

Define the following sequences 


f 9iU^ i ^ n 
\ P i = n. 


Thus for each i an H-measure Vi associated to vf equals Ofpu. 

On the other side we have that J^i '^r = (1 ~ (^n)xP -I- /", and the corresponding measure equals pu-\- Pf- 


Thus in order to constrain support of an H-measure by supports of corresponding components we have 
to impose additional assumptions on constituting sequences. More precisely, the following result holds. 
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Lemma 2.3. Let {9i) be an averaging sequence of positive numbers summing to 1, and let (u'^)n,i £ N 6e 
a family of uniformly bounded sequences, i.e. we assume there exists a constant Cu such that ||m”||l 2 < 

Cu,i, n £ N. 

Define the linear combination Vn = o-nd denote by fXi and v H-measures associated to (sub)se- 

quences (of) {u'(’)n and (vn), respectively. Then 

supp 1 /C Cl^ Ui supp/Xi^ . (2-10) 

Proof: Take an arbitrary pseudodifferential operator of order zero, P £ 4'°, with a symbol p(x, being 
compactly supported within the complement of the closure of U^supp fii. 

By the definition of H-measures we have 


« OO OO 

{n,p) =lim 


( 2 . 11 ) 


As P is a continuous operator on L^(R“) it follows that 


lim 


P 




OO OO 

6l,u”) (x) (x)d> 


< lim sup Cp 


L2 


< lim sup Cp j ^ ( 


( 2 . 12 ) 


where Cp is the bound of the operator P. The last sum is a remainder of a convergent series, and the 
above limit converges to zero uniformly with respect to n. 


Similarly, one shows the same property holds for lim„ 
we can exchange limits in (12.111) . getting 


/r. p{T.T (x) (Er Sfuj) (x)d> 


. Thus 


OO 


i=i j=i 




P {Oiu'f) (x) 6juf{x.)dK 


OO OO 

= ^ ^ ^ ,P) = 0, 

i=l j = l 


(2.13) 


where pij are H-measures determined by sequences (u") and [uf), supported within the closure given in 
(12.101) . outside which p is supported. □ 

As a consequence of the last lemma, in order to apply the localisation property within the analysis of 
observability estimates for solutions to (12.9L some kind of uniform boundedness on the solutions has to be 
assumed. 


Proposition 2.4. Suppose the observability inequality (ESI) holds for a solution ui to dsn. 

As for the system (ESI), suppose that norm of all the solutions Ui is dominated (up to a multi¬ 
plicative constant, independent of a choice of initial data) by the energy norm of ui. In addition as¬ 
sume that characteristic set {pi(t,x, r, ^) = 0} has no intersection with Cl^ Ui >2 {pi(f,x, r, £) = 0}^ for 

(t,x) £ (0,r) X w, (r, £ S'^, where pi stands for the principal symbol of the operator Pi. 

Then for any averaging sequence iOi) of positive numbers summing to 1, with 9i > 0, there exists a 
constant Cg such that the observability inequality 

Eii0)<Cgl [ [\J2(^^u^\^d^dt + \\(3o\\l-l + m\l -2 

\Jo Juj 

holds for any family of solutions {ui) to (I2S]) . 


(2.14) 
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Proof: Assume the contrary. Then there exist sequences of initial conditions (/3 q ), (/?"), and of associ¬ 
ated solutions (u"), such that 


1 = E^o) := mwh + > n 


/ 0 J uj 


+ 11 / 3^11 




(2.15) 


Let u be an H-measure associated to a (sub)sequence of Due to the inequality (12.151) . it equals 

zero on ( 0 , T) x oj. 

We split the last sum into two parts rewrite v in the form 

u = ui -b U2 + 25iui2, 


where ui and V 2 are H-measures associated to (sub)sequences (of) (0 ”m") and (X^^* respectively, while 

ui 2 is a measure corresponding to their product. In addition, ui = where by fii we denote a measure 

associated to a (sub)sequence (of) the i—th component u". 

From here the statement of the theorem is obtained easily (following the lines of the proof in finite 
discrete case. Proposition 12.ID . once we show that vi and U 2 have disjoint supports. 

By the localisation property for H-measures, each measure fii is supported within the set {pi(t, x, r, ^) = 

0 }. 

The assumption on the domination of solutions to (12.9D by an energy norm of ui, together with the 
constant initial energy £’"(0) implies uniform bound on solutions m”, both with respect to i and n. Thus we 
can apply Lemma 12.31 to conclude that 1/2 is supported within the set 

Cl( U >2 {Pz(t, X, T, I) = 0}) , 

which, due to the assumption on separation of the characteristics set, does not intersect the support of 
ui = 9lpi. As 9i is strictly positive, we get that it" converges to 0 strongly in L^((0,T) x uj), which 
contradicts the observability estimate ()2.2D . □ 


Remark 4. 

• The assumption of the last proposition requiring solutions Ui of ( 1 ^ to be dominated by the energy 

norm of ui occurs, for example, in a case of a system consisting of the operators of the same form. 
Pi = — Ai(t,x)^ • with uniformly bounded (both from below and above) coefficients and initial 

energies. 

The assumption on separation of characteristics sets in that case can be stated as 
Ai(t,x)|-| > (<) sup(inf) A,(t,x)| • I, (t, x) e (0, T) x w, ^ 7 b 0 , 

i >2 *>2 

i.e. the fastest (or the slowest) velocity is strictly separated from all the others. 

In that case the weak observability ()2.14|1 result is equivalent to the averaged controllability of the 
adjoint system up to a finite number of low frequencies. 

Of course, one can construct more general systems, including operators of different types as well, 
that satisfy the required boundedness assumption. 

• As already mentioned in previous subsection, obtaining corresponding strong observability result in 
this setting remains an open problem. 

• The constant Cg appearing in (I2.14|) can he taken uniformly for a family of averaging sequences, each 
satisfying 

(*) 9i>9,, 

(**) Er < efc, 

where 0 * € ( 0 , 1 ] and (sk) is a null sequence, both independent of a choice of a particular sequence 


2 . 3 . Simultaneous observability. The subsection deals with a problem of recovering energy of a system 
by observing an average of solutions on a suitable subdomain. For this purpose one has to estimate initial 
energies of all system components, unlike the case of the average observability where this was required just 
for the first one. A two component system is analysed firstly, while generalisations to a more dimensional 
case is discussed at the end. 
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We reconsider the system m assuming that P 2 is an evolution operator of the form 

P 2 = + A 2 , fceN, (2.16) 

where A 2 is an (uniformly) elliptic operator (in general different from Ai)^ and the problem for the pertur¬ 
bation U 2 is accompanied by a series of initial conditions 

{{dtyu2) (0) = e H-J(fl), 3=0,...,k-l. 

Its initial energy is denoted by 

fc-i 

£^2(0) = ^ ll 7 jllH-j(a) ■ 

0 

As in the previous subsection, we start with a weak observability result. 

Proposition 2.5. Suppose that there is a constant C, time T and an open subdomain uj such that for any 
choice of initial conditions the solutions to (EU) satisfy 

Eyo)<C [ [ \u^ydxdt, i = l,2. (2.17) 

Jo J u) 

In addition assume that characteristic sets {pi(t,x,r,^) = 0},* = 1,2 have no intersection for (t,x) S 
(0, T) X w, (r, ^) € where pi stands for the principal symbol of the operator Pi. 

Then for any 61,62 G R \ {0} there exists a constant Ce such that the observability inequality 

Ei{0) + £ 2 ( 0 ) < Cg ^ J l^iui-I--I-||/3o||jj-i + ||/3 i||jj-2 + ||7o|Ih-i + • ■ • + (2-18) 

holds for any pair of solutions {ui,U 2 ) to (12.11) . 

The result is obtained easily by following the steps of the proof presented above in the averaged ob¬ 
servability setting. Assuming the contrary and implying microlocal analysis tools, one shows that both 
components it" converge to 0 strongly on the observability region, thus obtaining the contradiction. 

However, a different approach is required in order to obtain the strong observability inequality for initial 
energy by removing compact terms in (12.1811 . It is based on a standard compactness-uniqueness procedure 
of reducing the observability for low frequencies to an elliptic unique continuation result mM- 

We introduce a subspace N{T) oi H = L^(H) x H“^(f2) x L^(H) x • • • x H^“^(r2), consisting of initial 
data for which the average of solutions to (HID vanishes on the observability region 

N{T) := {(/3o, /3i, 7 o,..., 7i-fe) S H\ 6 iUi + 62 U 2 = 0 on (0, T) x w}. 

Based on the relaxed observability inequality (I2.18P it follows that N{T) is a hnite dimensional space. 
Furthermore, the following characterisation holds. 

Lemma 2.6. We assume one of the following statements holds: 

a) The order k of time derivative in (12.1611 is odd. Coefficients of both the operators Pi and P 2 are time 
independent and of class C^’^, 

b) The time derivative order k is even, and — A 2 (or —(^J^^ — A 2 )) is an uniformly elliptic 
operator. Coefficients of both the operators Pi and P 2 are analytic. 

Then N{T) = {0}. 

Proof: One first shows that N{T) is an A-invariant, where A is an unbounded operator on H: 


/ 0 -1 

0 ^ 

Ai 0 


0 

1 

h- 1 

1 

to 

0 


with the domain D{A) = Hj(fl) x L^(fl) x H(,(H) x • • • x ^(^)- 


(2.19) 
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Being ^-invariant and finite-dimensional, it contains an eigenfunction of A. Thus there is a A £ C and 
(, 80 ,^ 1 , 70 , ■ ■ -,71-^) G such that 


^i/3o — 

AI270 = (—1)^ ^-^^70 

/3i = — A/3o 


( 2 . 20 ) 


7i = (-!)■’^-^70, j = 


By the definition of N{T) it follows OiUi + O 2 U 2 = 0 on (0,T) x w, and specially 

+ ^27 o = 0 on w . (2.21) 


At this level, we want to show that each assumption of the lemma implies /3o = 7o = 0. 

a) As Ai, i = 1, 2 are positive operators and k is odd, from (12.201) it follows that one of functions Po, 70 is 
trivial. By relation (12.211) it follows that the other one also equals zero on ui. Being an eigenfunction 
of an elliptic operator, the unique continuation argument (e.g. [5J Theorem 3]) implies it is zero 
everywhere. 

b) Analyticity of coefficients implies analyticity of eigenfunctions. Specially it follows 0i/3o + ^270 = 0 
everywhere, and relations (12.201) imply 

(AJ/' - A2)Po = 0. 

k /2 

Assumptions on the operator Ap — A 2 imply /3o = 0 on fl. 


□ 


Remark 5. In a special case A 2 = —div(c2(x)V) one easily proves that the last Lemma holds with analytic 
coefficients ci,C 2 being separated just on an arbitrary non-empty open set, and not on the entire LI. 

Theorem 2.7. Under the assumptions of Provosition \ 2.1\ and Lemma \2.6\ there is a positive constant Cg 
such that the strong observability ineguality holds: 

Ei{0) + E2{0 )<Cg [ [ \9iUi+e2 U2\^dxdt. (2.22) 

Jo J uj 

Proof: As in the proof of Proposition 12.11 let us suppose the contrary. Then there exists a sequence of 
solutions uf,U 2 to (12.11) such that if”(0) = 1 and 

\9iu'f -f 92U2\'^d'Kdt —0. 

Thus for weak limits (mi,M 2 ) of solutions on the observability region we have 9iUi 92 U 2 = 0, implying 
(mi(0), 9t'Ui(0), 1 / 2 ( 0 ),... , {{dt)^~^U 2 ) (0)) £ N(T). By means of the above lemma and taking into account 
the relaxed observability inequality, it follows 

1 < Cg ^ + 92 U^\^dKdt + ||/3o"||^-i + + I|7o”IIh- + • ■ • + Il7fc”-lll|-.^ 0, 

thus obtaining a contradiction. □ 

We close this subsection by the following remarks. 

Remark 6. 

• If P 2 is a second order evolution operator the strong observability ineguality (12.221) is equivalent to 
the simultaneous controllability of the adjoint system, also studied in [10) . by which one controls each 
component individually (and not just their average). 

• The notion of simultaneous observability is stronger than the average one, as it estimates energy of 
all system components, whose initial data, in this case, are not related. Consequently, it requires 
stronger assumption of GCC being satisfied by each component. 
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• The application of the compactness-uniqueness procedure in the passage from the weak to the strong 
observability estimate allows the perturbation P 2 to be an evolution operator with an arbitrary elliptic 
part. However, such approach is not possible in the averaged observability setting. Namely, in order 
for subspace N{T) to be finite dimensional one has to relate the initial data of two components by a 
bounded linear operator. But such constrain would not be preserved under action of the operator A 
given by (j2.19|) . and as a consequence N(T) would not be A-invariant. 

• The weak observability result 112.181) is easily generalised to a system with a finite number of compo¬ 
nents, under assumption that the characteristic sets of all operators are mutually disjoint. 

• As in the averaged observability case, the generalisation of strong estimate result (j2.22|l to a system 
consisting of more than two components has still not been obtained, and is a subject of the current 
investigations. 

3. Observation of the Schr5dinger equation under uncertain perturbations 

In this section we consider a system in which the first component, the one whose energy is observed, 
satisfies the Schrodinger equation, while the second one, corresponding to a perturbation, is governed by an 
evolution operator P 2 : 

Piui = idtui + div (Ai(x)Vui) = 0, (i,x) € x fl 

P 2 U 2 = 0, {t, x) € R+ X n 

(3-1) 

ui = 0, (i,x) s R'*' X dfl 

uii 0 ,-) = 

As in the study of perturbations of the wave dynamics in Section [2l we specify no initial or boundary 
conditions for the second operators, we just assume that corresponding problem is well posed and that it 
admits an solution. As for the system coefficients, as before we impose merely boundedness and continuous 
assumptions, and suppose that Ai is a positive definite matrix field. 

3.1. Averaged observability under non-parabolic perturbations. For the reasons explained below, 
in this subsection we restrict the analysis to evolution operators P 2 of order strictly larger than one. In that 
case the stability of the Schrodinger observability estimate is given by the next theorem. 

Theorem 3.1. Suppose that there is a constant C, time T and an open subdomain lo such that for any 
choice of initial datum the solution ui of (EH) satisfies 

Ei{0)-.= 1^01112 <C [ [ \ui\^d^dt. (3.2) 

Jq J u) 

In addition, for the system EH we assume the following: 

a) The perturbation operator is an evolution operator of the form (12.51) and of the order k > 1. 

b) The initial values of solutions Ui,i = 1,2 are related by a linear operator such that whenever 
((6IiMi(0) + 02 U 2 iO))\u = 0) then (ui(0)|<^ = M2(0 )||c,; = O). 

Then, for any 61 G (0,1] there exists a constant Ce such that the observability inequality 

Ei{Q)<Cg f f \9iUi-\-92U2\^dxdt (3.3) 

Jo J u) 

holds for any pair of solutions (ui,U 2 ) to (12.11) . 

The proof goes similarly as for the observations of the wave equation. Required conditions a), b) are 
necessary for obtaining the strong observability inequality, without a compact term. 

On the other hand, in order to obtain a relaxed inequality with a compact term, no assumption is 
required at all. Namely, the assumption on separation of characteristic sets {pi(t,x, r, ^) = 0} required 
in Proposition 12.II becomes superfluous in this setting, as being directly satisfied by an arbitrary evolution 
operator P 2 of order k strictly larger than 1. Namely, its characteristic set does not contain the poles 4 = 0 
of the unit sphere in the dual space, which constitute the characteristic set of the Schrodinger operator Pi. 
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However, every Schrodinger or the heat operator P 2 fails to satisfy the assumption. Namely, no matter 
the coefficients entering the equation, both have characteristic set within the poles ^ = 0, same as Pi, and 
localisation principle fail to distinguish corresponding H-measures. To analyse such a system one needs 
a microlocal tool better adopted to a study of parabolic problems. Namely, original H-measures were 
constructed with the aim of analysing hyperbolic problems and are not capable to distinguish differences 
between the time and space variables that are intrinsic to parabolic equations. Their variant, parabolic 
H-measures, was recently introduced with the purpose of overcoming the mentioned constraint. 

3.2. Parabolic H-measures. Parabolic H-measures were first introduced in [I], while more exhausted 
introduction can be found in [2] , elaborating in particular their basic properties: localisation and propagation 
one. The former will be used in the next subsection for proving the stability of the observability estimate 
for the Schrodinger operator under parabolic perturbation. 

Here we expose the basic result on parabolic H-measures used in the note. 

The main idea in their construction is to replace the projection along the straight rays in the dual space, 
determined by the term 4/|4| in the definition of the original H-measures, by the one going along meridians 
of paraboloids t = a$, . The hypersurface on which the dual space (except the origin) is projected is a 
rotational ellipsoid 



Although the ellipsoid might seem as an unnatural choice of the surface on which one want to construct 
the parabolic H-measures, crucial in its choice was that the curves along which the projections are taken 
intersect it in the normal direction, as it was in the classical case, where the rays radiating from origin are 
perpendicular to the unit sphere. The mentioned normality property enables study of propagation properties 
of the measures. 

Concerning the applications, the most important is that the new tool is also a kind of defect measure, 
in the sense that null parabolic H-measure is equivalent to strong convergence of (v„) in 

In order to formulate the localisation principle, we first introduce some special anisotropic (Sobolev) 
function spaces 

Hi.*(Ri-Ld) .= g 5' : g L2(R1+'^)} ,s g R, 

where fcp(T,^) := + (27rr)^ -|- (27r|4|)‘^ is the weight function. These are Hilbert spaces and they are 

particular examples of more general Hormander spaces Bp^k described in [9j §10.1]. 

We also define the fractional derivative: \fdt as a pseudodifferential operator with a polyhomogeneous 
symbol i.e. 

\fdtu = T ^•\/27rzTu(T)^ . 

Theorem 3.2. (localisation principle for parabolic H-measures, [2]) Let (v„) be a sequence of func¬ 
tions uniformly compactly supported in t and converging weakly to zero in L^(R^“'"‘^; C’'), and let for s S N 

V^'(AV)+ ^ a“(A“v„)^0 strongly m , 

\a\=s 

where A*^, A“ are continuous and bounded matrix coefficients, while a G Nq. 

Then for the associated parabolic H-measure p we have 



In particular, the principle implies that the measure p is supported within the parabolic characteristic 
det {[yf 2'n:iTyAf) + ^ (27ri4)“AQ,^ = 0, (t, |) e P 

|a|=s 


set. 
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Example 2. (Application of the localisation principle to various equations) 

In all the examples we assume that the equation coefficients satisfy the assumptions of the preceding 
theorem, i.e are continuous and bounded. 

• The Schrodinger equation 

Let (u„) be a sequence of solutions to the Schrodinger equation 

idtUn + div (A(t, x) Vm„) = 0. 

If (un) is bounded in L2(R+ X R'^) then the associated parabolic H-measure p. satisfies 

{ 2 'kt + 4?!^A(t, x)| • = 0, 

implying p is supported in points of the form 2 t:t = —4?!^ A(t, x)^ • 

• The heat equation 

Let (un) be a sequence of solutions to the heat equation 

dtUn - div (A(t, x)Vu„) = 0 

where A is a bounded, positive definite matrix field. If {un) is bounded in L^(R''' X R"^) then the 
associated parabolic H-measure p satisfies 

(27riT + 471^A(t, x)^ • ^)p = 0. 

As (parabolic) H-measures live on a hypersurface in the dual space excluding the origin, and A is 
positive definite, the term in the braces above never equals zero, implying p is a trivial (null) measure. 

• The wave equation 

Let (un) be a sequence of solutions to the wave equation 

dttUn - div {A{t, x)Wun) = 0. 

If {un) is bounded in L2(R+ X R"*) then the associated parabolic H-measure p satisfies 

Att'^t^P = 0 , 

implying p is supported in the equator t = 0 of the hypersurface P. 


From the given examples it is clear that by taking two Schrodinger operators with separated coefficients 
we shall be able to distinguish corresponding parabolic H-measures. Of course, the distinction is also possible 
if the considered operators are of different type (e.g. the Schrodinger and the wave one), as it was the case with 
the original H-measures. This enables the generalisation of Theorem 13.II represented in the next subsection. 


3.3. Averaged observability under parabolic perturbations. We reconsider the system (13.11) with the 
aim of obtaining the stability of the observability estimate for the Schrodinger operator under perturbation 
determined by the operator P 2 under minimal conditions on the latter. In particular, we want allow it to be 
the Schrodinger, or some other parabolic type operator, like Pi. 

The results of previous subsection, in particular the localisation property (Theorem 13.2p provide the 
following result. 


Theorem 3.3. In addition to assumptions of Theorem \.S.l\ we allow P 2 to be an evolution operator of order 
one, and assume the separation of coefficients C 2 (x) ^ —i holds on uj in the case k = 1. 

Then for any 0i € (0,1] there exists a constant Cg such that the strong observability inequality 


Ei{Q)<Ce f [ \9iuiP 92U2\^dJ^dt 

Jo J u) 


(3.4) 


holds for any pair of solutions {ui,U 2 ) to ()3.1I1 . 


The proof goes along the same lines as the one of the Theorem 12.21 using that the parabolic H-measures 
share the basic properties (positive definiteness, diagonal domination) with their original (hyperbolic) coun¬ 
terparts. Crucial in the proof is separation of corresponding parabolic characteristic sets. For an evolu¬ 
tion operator P 2 of order k > 2 that set is restricted to the equator r = 0 of the hypersurface P, where 
Pi = 27rT -I- 47r^A(t,x)^ • ^ differs from zero (cf. Example [2]). For k = 1 the separation follows by the 
assumption C2(x) 7^ —i. 
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The novelty obtained by the application of parabolic measures is that both the operators entering the 
system are allowed to be of the same type (e.g. two Schrodinger operators with coefficients being 

separated on the observability region (0,T) x uj). In that case the observability estimate (13.41) is equivalent 
to the averaged control of the adjoint system. 

Of course, operators of different types are admissible as well, thus the last theorem incorporates the 
results of Subsection [3T] as well. 

Finally, let us mention that as in Section [2] one can obtain analogous result for the robust observability 
of the Schrodinger equation in the simultaneous and infinite discrete setting. 

4. Conclusion 

One of the most interesting aspects of the above results refers to a system determined by two operators 
of a same type - either the Schrodinger or the wave one. The weak observability estimates in that case 
correspond to the average control of solutions to adjoint system up to a finite number of low frequencies. The 
required assumptions are optimised: the coefficients are merely continuous and separated on the observability 
region. The strong averaged observability result, corresponding to the exact averaged controllability, assumes 
in addition initial data of two components to be related in an appropriate manner. Its proof, unlike in 
m, employs neither propagation property of H-measures nor unique continuation procedure, therefore not 
requiring additional smoothness assumptions. 

The paper restricts to decoupled systems, but one could analyse more general ones with coupling in 
lower order terms (e.g. m)- Such terms do not effect microlocal properties of solutions, thus enabling a 
generalisation of the obtained estimates. 

Surely of interest would be a generalisation of the results to a parameter dependent system with a para¬ 
meter ranging over a continuous set. A problem has been analysed in m, exploring the heat and Schrodinger 
equations with a randomly distributed parameter. The constant coefficients operators are considered, with 
eigenfunctions of corresponding elliptic parts being independent of the parameter value. Thus the explored 
averages are presented as a solution (or a superposition of two solutions) to a similar evolution problem(s), 
which is crucial in the proof. 

As the next step in that direction, one would consider a system of equations determined by evolution 
operators whose elliptic parts coincide up to a scalar function: 

P(u) = -t c(x, u)A , 

with u being a parameter, and A the elliptic part of P. The variable dependence of the coefficient c would 
require another techniques from those applied in m- In this paper we have obtained a corresponding result 
for a parameter ranging over an infinite discrete set, but only at the level of a weak observability estimate 
with compact terms. 

By using transmutation techniques developed in [^, the simultaneous observability result of subsection 
12.31 can be employed in order to obtain controllability and observability properties for a system of heat 
equations (cf. [10)1. Similar transmutation procedure can be constructed, transferring Schrodinger into wave 
type problems, and vice versa. Its application will result in observability estimates for a system of Schrodinger 
equations, derived from the corresponding results for the wave ones. It would be interesting to compare such 
obtained results with those obtained directly in Section[3]by means of parabolic H-measures, and to compare 
the efficiency of methods applied by each approach. 
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